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Abstract 



A search for the rare decays Bg — > }i }i^ and B — > }i }i^ is performed in pp 
collisions at ^s = 7 TeV, with a data sample corresponding to an integrated lu- 
minosity of 1.14 fb^ , collected by the CMS experiment at the LHC. In both cases, 
the number of events observed after all selection requirements is consistent with 
expectations from background and standard-model signal predictions. The result- 
ing upper limits on the branching fractions are B(B^ — ?> ji^jx^) < 1.9 x 10^^ and 
B(B° -^ }i^H~) < 4.6 X 10"^ at 95% confidence level. 
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In the standard model (SM) of particle physics, flavor-changing neutral current (FCNC) decays 
are forbidden at tree level and can only proceed through higher-order loop diagrams. The 
decays B(i(s) — ?► £+£^ (where I = e, }i), besides involving b — ?► s(d) FCNC transitions through 
Penguin and box diagrams, are helicity suppressed by factors of {nii/m-of-, where m( and wb 
are the masses of the lepton and B meson, respectively. They also require an internal quark 
annihilation within the B meson that further reduces the decay rate by (/b/wjb)^, where /b is 
the decay constant of the B meson. 

The SM-predicted branching fractions, 5(B° -^ }i^}i') = (3.2 ± 0.2) x 10"'' and B{B° -^ 
}i^fi^) = (1.0 ± 0.1) X 10^^*^ [IJ, are significantly enhanced in several extensions of the SM, 
although in some cases the decay rates are lowered ||2l. For example, in the minimal supersym- 
metric extension of the SM, the rates are strongly enhanced at large values of tan/5 ||3[|11. In 
specific models involving leptoquarks fS] and in supersymmetric models with non-universal 
Higgs masses [6|, the Bg — > fi^ji^ and B*^ — > ji^ji^ branching fractions can be enhanced by 
different factors and, therefore, both channels must be studied in parallel. Several experiments 
have published upper limits at 95% confidence level (CL) on these decays: S(Bg -^ H^h) < 
5.1 X 10-^ by DO [Zi; B{B° -^ ^^^') < 5.8 x 10"^ and B(B^ -^ H^V') < 1-8 x 10"^ by 
CDF [8J; i3(B0 -^ ^+j^-) < 5.6 x 10"^ and B{B° -^ ^+^-) < 1.5 x 10"^ by LHCb l9|. CDF 
recently reported a new limit of B{B^ -^ F^H) < 6.0 x 10^^ and an excess of Bg — > ji^pi^ 
events, corresponding to 6(B0 -^ /^+f/") = (1-8+^:9) ^ ^^^^ ESI- 

In this Letter, a simultaneous search for the Bg — > }i^}i^ and B^ — > ji^ji^ decays is presented, 
using a data sample of pp collisions at \/s = 7 TeV, corresponding to an integrated luminosity 
of (1.14 ± 0.07) fb"^ collected in the first half of 2011 by the Compact Muon Solenoid (CMS) 
experiment at the Large Hadron Collider (LHC). An event-counting experiment is performed 
in dimuon mass regions around the Bg and B*^ masses. To avoid any possible bias, the sig- 
nal region w^as kept blind until after all selection criteria were established. The backgrounds 
are evaluated from the yields measured in data mass sidebands and from Monte Carlo (MC) 
simulations for rare hadronic two-body B decays. The MC event samples are generated with 
Pythia 6.409 ini, the unstable particles are decayed via EvtGen [12], and the detector re- 
sponse is simulated with Geant4 [13]. Events of the type B^ —> ]/ip¥i^, ]/xp -^ F^V-^ are used 
as a normalization sample to minimize uncertainties related to the bb production cross section 
and to the integrated luminosity. The signal and normalization efficiencies are determined with 
MC simulation studies. A control sample of reconstructed B*^ — )■ ]/'ip(p, ]/tp -^ }i^}i^ events is 
used to validate the MC simulation (such as the Bg transverse momentum px spectrum) and to 
evaluate potential effects resulting from differences in fragmentation between B+ and Bg. The 
analysis is not affected by multiple pp collisions in the same bunch crossing (pileup) because 
the spatial vertex resolution is good enough to correctly identify the pp vertex from which 
signal candidates originate. In the present data set, an average of 5.5 primary vertices are re- 
constructed per event. 

A detailed description of the CMS experiment can be found in Ref. [^4^. The main subdetectors 
used in this analysis are the silicon tracker, composed of pixel and strip layers immersed in a 
3.8 T axial magnetic field, and the muon stations, made of gas-ionization detectors embedded 
in the steel return yoke, and divided into a barrel section and two endcaps. The muons are 
tracked w^ithin the pseudorapidity region |f/ 1 < 2.4, where ?/ = — ln[tan(0/2)] and d is the polar 
angle with respect to the counterclockwise beam direction. A muon px resolution of about 
1.5% is obtained for muons in this analysis. The events are selected w^ith a two-level trigger 
system. The first level only requires two muon candidates, without an explicit px requirement, 
while the high-level trigger (HLT) uses additional information from the silicon tracker. The 
HLT selection for the signal data sample requires two muons each with px > 2 GeV, dimuon 



Pt > 4 GeV, invariant mass within 4.8 < m^^ < 6.0 GeV, and a three-dimensional (3D) distance 
of closest approach to each other d'^^ < 5 mm. 

The normalization (B^ — ?► J/t/^K^) and control (B° — )■ ]/^(p) samples were collected w^ith HLT 
requirements gradually tightened as the LHC luminosity increased. This time evolution does 
not affect the analysis presented here, which uses selection criteria significantly tighter than 
any trigger requirements. More than 95% of the normalization and control sample events w^ere 
collected by requiring two muons each with px > 3GeV, dimuon pj > 6.9 GeV, invariant 
mass within 2.9 < m^,^ < 3.3 GeV, d'^^ < 5 mm, and a larger than 0.5% probability of the x^ 
per degree of freedom (dof) of the dimuon vertex fit. Two additional trigger requirements, 
measured in the transverse plane, significantly reduce the rate of prompt ]/tp candidates: the 
significance of the flight distance ixy/^i^xy) rnust be larger than 3, where £xij is the distance 
between the primary and dimuon vertices and cr[ixy) is its uncertainty; and the pointing angle 
ocxy between the B candidate momentum and the vector from the primary vertex to the dimuon 
vertex must fulfill cos ocxy > 0.9. The average trigger efficiency for events in the signal and 
normalization samples is about 80%, as determined from MC simulation. The uncertainty on 
the ratio of trigger efficiencies between the signal and normalization samples is estimated to be 
2% by comparing these ratios in simulation studies and in data. 

Muon candidates are required to be reconstructed by two different algorithms, one matching 
silicon-tracker tracks to segments in the muon stations, and the other performing global fits 
using tracks in both detector systems IflSl . The uncertainty on the ratio of muon identification 
efficiencies between the signal and normalization samples is estimated to be 5%. 

The B — ?► ji^fi^ candidates require two oppositely charged muons with an invariant mass 
in the region 4.9 < m^^, < 5.9 GeV, after constraining their tracks to come from a common 
vertex. The B candidate momentum and vertex position are used to choose a primary vertex 
based on the distance of closest approach. Since the background level depends significantly 
on the pseudorapidity of the B candidate, the events are separated into two categories: the 
"barrel channel" contains the candidates where both muons have |7/| < 1.4 and the "endcap 
channel" contains those where at least one muon has \ri\ > 1.4. An isolation variable I = 
Pt(B) / (Pt(B) + Litrk Pt) is calculated from the transverse momentum of the B candidate Pt{B) 
and the transverse momenta of all other charged tracks satisfying AR = \/{ArjY + {^<pY < 1/ 
where A;/ and A(p are the differences in pseudorapidity and azimuthal angle between a charged 
track and the B candidate momentum. The sum includes all tracks with px > 0-9 GeV that are 
consistent with originating from the same primary vertex as the B candidate or have a distance 
of closest approach dca < 0.5 mm with respect to the B vertex. The minimum distance of closest 
approach with respect to the B vertex among all tracks in the event, d^, is also determined 
as a complementary isolation variable. Figure IT] illustrates the transverse momentum, the 3D 
pointing angle ^sd, the 3D flight length significance ^sd/c^I^sd)/ and the isolation distributions 
for signal MC and for sideband background data events. The sideband covers the range 4.9 < 
m^ip < 5.9 GeV, excluding the signal window 5.2 < m^,^ < 5.45 GeV. 

The following selection requirements were optimized for the best expected upper limit using 
MC signal events and data sideband events. The requirements were established before ob- 
serving the number of data events in the signal region. The optimized requirements include 
Pt > 4.5 GeV on one muon and px > 4.0 GeV on the other, B candidate px > 6.5 GeV, I > 0.75, 
and B-vertex fit x^/doi < 1.6. Two requirements are different for the barrel and endcap chan- 
nels: asD < 0.050 (0.025) and £3D/o"(4d) > 15.0 (20.0) for the barrel (endcap). Furthermore, 
for events in the endcap there is an additional requirement, d^™ > 0.15 mm. The signal effi- 
ciencies £tot of these selections are provided in Table [l] The dimuon mass resolution for signal 



Table 1: The event selection efficiencies for signal events Stot, the SM-predicted number of signal 
events N^^^j, the expected number of combinatorial background events N^^^]^ and peaking 

background events N^^^, and the number of observed events Nobs in the barrel and endcap 

//+/^^ and B*^ — > }i^H^ ■ 



channels for Bg 





Barrel 


Endcap 




bo -^ ^+^- 


bo -^ n+ji- 


bo -^ ^+^- 


B° ^ }i+V- 


£tot 


(3.6 ±0.4) xlO-3 
0.065 ±0.011 


(3.6 ±0.4) xlO-3 
0.80 ±0.16 


(2.1 ±0.2) X 10-3 
0.025 ± 0.004 


(2.1 ±0.2) xlO-3 
0.36 ± 0.07 


''^peak 


0.40 ± 0.23 
0.25 ± 0.06 


0.60 ± 0.35 
0.07 ±0.02 


0.53 ± 0.27 
0.16 ±0.04 


0.80 ± 0.40 
0.04 ± 0.01 


Nobs 





2 


1 


1 



events depends on the pseudorapidity of the B candidate and ranges from 36 MeV for tj 
85 MeV for \rj\ > 1.8, as determined from simulated signal. 

The reconstruction of B± -^ ]/tpK^ ^ m+u"K± (B° 



0,to 



]/tp(p -^ }i^}i K+K ) candidates re- 



quires two oppositely-charged muons with an invariant mass in the range 3.0-3.2 GeV, which 
are combined with one (two) track(s), assumed to be (a) kaon(s), fulfilling pj > 0.5 GeV and 
\r]\ < 2.4. To ensure a well-measured trigger efficiency, the selected candidates must have 
dimuon pj > 7 GeV and the two muons must bend away from each other in the magnetic field 
(to avoid spurious detector-induced pair correlations). The d^^ between all pairs among the 
three (four) tracks is required to be less than 1 mm. For B*^ — ?► ]/xp(p candidates the two assumed 
kaon tracks must have an invariant mass in the range 0.995-1.045 GeV and AR(K+, K^ ) < 0.25. 
The tracks from all decay products are used in the B-vertex fit and only B candidates with an 
invariant mass in the range 4.8-6.0 GeV are considered. The efficiencies of individual selection 
criteria agree to better than 4% (6%) between data and MC simulation for the normalization 
(control) sample. Figure compares several distributions for B" — ?► ]/tp(p candidates between 
MC simulation and sideband-subtracted data. 



-> ]A^jA K+, including the detector acceptance, is efot = 
^, respectively for the barrel and endcap channels, where 



The total efficiency for B^ -> ]/ipK^ 
{7.7 ± 0.8) X 10"^ and (2.7 ± 0.3) x 10 
statistical and systematic uncertainties are combined. The acceptance has a systematic uncer- 
tainty of 4%, estimated by comparing the values obtained with different bb production mecha- 
nisms (gluon splitting, flavor excitation, and flavor creation). The uncertainty on the event se- 
lection efficiency for the B^ — > J/t/^K^ normalization sample is 4%, evaluated from differences 
between measured and simulated B^ — > J/t/^K^ events. The uncertainty on the signal efficiency 
(7.9%) is evaluated using the B*^ — > ]/ip(p control sample. The invariant mass distributions are 
fitted w^ith a Gaussian function for the signal and an exponential (barrel) or a first-degree poly- 
nomial (endcap) plus an error function for the background, as shown in Fig. [3] Applying the 
same selection requirements as for the signal sample, the observed number of B^ — > J/t/^K^ 
candidates in the barrel (endcap) channel is N^^^ = 13045 ± 652 (4450 ± 222). The uncertainty 
includes a systematic term caused by fit and background parametrization effects, estimated to 
be 5% from MC studies. 

To quantify a possible dependence on the pileup, the efficiencies of the isolation and the flight 
length significance requirements are calculated as functions of the number of reconstructed 
primary vertices. No dependence is observed for events with up to 12 primary vertices for the 
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Figure 1: Comparison of MC signal and sideband data distributions, for the transverse mo- 
mentum (top left), the 3D pointing angle (top right), the flight length significance (bottom left), 
and the isolation (bottom right). The MC histograms are normalized to the number of events 
in the data. 



CMS, 1.14 fb^^ 



s = 7 TeV 




CMS, 1.14 fb 



s = 7 TeV 



CMS, 1.14 fb"^ 



30 40 50 

P^(B) [GeV] 

Vs = 7 TeV 



Data 

B° ^ J/>(/ 4) (MC) 




■§2000 



CO 

O 



1500 



1000 



500 




CMS, 1.14 fb"^ 



s = 7 TeV 




80 100 0.2 0.4 0.6 0.8 1 

'3d/*^('3d) isolation 

Figure 2: Comparison of measured and simulated B° — )■ J/)/?(/' distributions, for the transverse 
momentum (top left), the 3D pointing angle (top right), the flight length significance (bottom 
left), and the isolation (bottom right). The MC histograms are normalized to the number of 
events in the data. 



normalization and control samples. 



The Bg — ?► |^^^ branching fraction is measured separately in the barrel and endcap channels 
using 



6(Bl^ 



H ^ 



Ns fuc 



N^. /^ 



s £tot 



(1) 



and analogously for the B*^ — ?► ji^fi case, where Ns is the background-subtracted number of 
observed B(j(s) — )■ }i~^}i~ candidates in the signal window (5.3 < m^,, < 5.45 GeV for Bg and 
5.2 < nifiif < 5.3 GeV for B°) and etot is the total signal efficiency of all selection requirements. 
The ratio of the Bg and B+ meson production fractions is fs/fu = 0.282 ± 0.037 and 6(B+) = 
B{B+ -^ ]/xpK+ -^ ^+|^-K+) = (6.0 ± 0.2) x 10"^ 116J. (We use /s = 0.113 ± 0.013 and 
/u = 0.401 ± 0.013 from the main section of Ref. ||T6| and account for the correlations in the 
ratio.) 

Events in the signal window can result from real signal decays, combinatorial background, and 
"peaking" background from decays of the type B^tg\ — ?► hh', where h, h' are charged hadrons 
misidentified as muons. The expected number of signal events, N^g^ai' ^^ calculated assuming 
the SM branching fraction and is normalized to the B+ yield. The expected number of combi- 
natorial background events, N^^^^ = 3(4) in the barrel (endcap), is evaluated by interpolating 
to the signal window the number of events observed in the sideband regions. The interpolation 
procedure assumes a flat background shape and has a systematic uncertainty of 4%, evaluated 
by varying the flight length significance selections and by using a floating slope. The expected 
number of peaking background events, N^^^j,, is evaluated from MC simulation and muon 

misidentification rates measured in Kg — >• 71+ tt^, (p — >• K+K~, and A — )■ pw samples [15|. 
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Figure 3: B+ — ?► J/i/^K+ invariant mass distributions in the barrel (left) and endcap (right) chan- 
nels. The solid (dashed) lines show the fits to the data (background). 

Figured shows the measured dimuon invariant mass distributions. Three events are observed 
in the B^ — > }i^}i^ signal windows (two in the barrel and one in the endcap), while only one 
event is observed in the B^ — > }i^}i^ endcap channel. This observation is consistent with the 
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Figure 4: Dimuon invariant mass distributions in the barrel (left) and endcap (right) channels. 
The signal windows for Bg and B*^ are indicated by horizontal lines. 

SM expectation for signal plus background. Upper limits are determined with the CLg ap- 
proach HZl. Table [T] shows the values needed for the extraction of the results, separately for 
the barrel and endcap channels. The obtained upper limits on the branching fractions are 
5(b0 -^ /^+^-) < 1.9 X 10-^ (1.6 X 10-8) and 6(B° -^ f/+/^-) < 4.6 x 10"'^ (3.7 x 10"''), 
at 95% (90%) CL. The median expected upper limits at 95% CL are 1.8 x 10"^ (4.8 x 10"'^) 
for Bg -^ H^j^^ (B" -^ ji^pi^). The background-only p value is 0.11 (0.40) for Bg -^ j^^H^ 
(B^ — > ji^ji^), corresponding to 1.2 (0.27) standard deviations. The p value is 0.053 when as- 
suming a Bg — >• }i^}i~ signal at 5.6 times the SM value, as reported in Ref. IITOl . 



In summary, a search for the rare decays Bg 



¥^¥ 



andBO 



^+}l- 



has been performed 



on a data sample of pp collisions at yjs = 7 TeV corresponding to an integrated luminosity of 
1.14 fb^ . The observed event yields are consistent with those expected adding background 
and SM signals. Upper limits on the branching fractions have been determined at 90% and 

95% CL. 
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